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ABSTRACT 

X-ray observations are a powerful diagnostic tool for transport, acceleration, 
and heating of electrons in solar flares. Height and size measurements of X- 
ray footpoints sources can be used to determine the chromospheric density and 
constrain the parameters of magnetic field convergence and electron pitch-angle 
evolution. We investigate the influence of the chromospheric density, magnetic 
mirroring and coUisional pitch-angle scattering on the size of X-ray sources. The 
time-independent Fokker-Planck equation for electron transport is solved numer- 
ically and analytically to find the electron distribution as a function of height 
above the photosphere. From this distribution, the expected X-ray flux as a 
function of height, its peak height and full width at half maximum are calculated 
and compared with RHESSI observations. A purely instrumental explanation 
for the observed source size was ruled out by using simulated RHESSI images. 
We find that magnetic mirroring and coUisional pitch-angle scattering tend to 
change the electron flux such that electrons are stopped higher in the atmo- 
sphere compared with the simple case with coUisional energy loss only. However, 
the resulting X-ray flux is dominated by the density structure in the chromo- 
sphere and only marginal increases in source width are found. Very high loop 
densities (> lO^^cm"^) could explain the observed sizes at higher energies, but 
are unrealistic and would result in no footpoint emission below about 40 keV, 
contrary to observations. We conclude that within a monolithic density model 
the vertical sizes are given mostly by the density scale-height and are predicted 
smaller than the RHESSI results show. 



Subject headings: Sun: flares - Sun: X-rays, 7-rays - Sun: Chromosphere 
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Introduction 



The transport of flare-accelerated electrons and the generation of hard X-ray (HXR) 
emission in the solar atmosphere is one of the most important and widely used diagnostics of 
flare accelerated electrons. These supra-thermal particles precipitate along the field lines of 
a magnetic loop from the acceleration site towards the denser regions of the chromosphere. 
They undergo Coulomb collisions with electrons and ions in the ambient plasma, and can 
be observed via their bremsstrahlung emission. The bulk of the observed HXR emission 
comes from the footpoints of magnetic structures where the density is high and electrons 
lose their energy completely, with the electron stopping location (depth) determined by the 
initial electron energy and the ambient density. 



Using this energy dependency in combination with observations from RHESSI (ILin et al. 



2OO2I ). it has become possible to characterize the structure of X-ra y sources. Assuming col- 



lisional transport, the chromospheric density could be inferred (Aschwanden et al 



Kontar et al 



2008L I2OIOI : ISaint-Hilaire et al.ll2010l : iBattadia fc Kontar 



2002 



2OIII). These obser- 



vations suggest that the bulk of HXRs in the range 30-100 keV is produced at heights of 
~ 0.7 — 1.2 Mm. The density scal e heights generally a gree well with hydrostatic chromo- 
spheric den sity models such as (e.g. IVernazza et al.lll981l) . However, using RHESSI visibility 



techniques (IHurford et al.l 12002 



Schmahl et al. 



20071), it is also possible to infer the energy 



deperi d ent source sizes with bette r than arcsecond accuracy as shown by iBattaglia fc Kontar 
( 120111 ): iKontar et al.l ( l2008l . |2010( ). The source sizes in the direction aloi ig the magnetic field 
inferr ed by these authors are in the range from 2 arcsec to 6 arcsec ( IBattaglia fc Kontar 
201 ll ) which is up to a factor of 4 larger than what would be expected frorn coUisi onal trans- 
port in the same density profile. To explain this discrepancy, iKontar et al.l (120 lOl ) suggested 
that the electrons propagate in a multi-thread loop with different density profiles along each 
thread, so that the X-ray source positions are the same as in the case of collisional transport 
in a single densi ty loop, but the ver t ical s ize is enlarged. However, even within a single 
monolithic loop, iBattaglia fc Kontarl ( 1201 ll ) indicated a number of processes which might 
increase the vertical extent of HXR sources. 

To investigate those effects quantitatively one has to solve the equations for electron 
transport in the footpoints. The electron transport problem has been considered in the past 
by several authors to investigate various aspects of electron transport, electron trapping and 
energy losses in the solar atmosphere. Both s e mi-an a lytical and full nuni e rical s oluti o ns ca n 
be found in the lite r ature . iLeach &: PetrosianI (Il981l): Kovalev fc KorolevI ( 1l98ll ) ; iBail ( 1l982l ) ; 



McClementd (Il990l . |l992|), and 



Kocharov et al.l (I2OOOI ) use numerical methods to investigate 



the ef fect of the magnetic fi eld g eometry and magnetic trapping on the electron spectrum 
while iMacKinnon fc Craigl (Il99ll ) use a test-particle method of solution, applying it to a 
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situation without magnetic mirroring. Many of those studies focused on the electron trans- 
port and the effect of scattering and magnetic mirroring on the electron distribution. To 
compare these models with actual obser vations, one has to go a st ep further and model the 
resulting X-ray flux, as was done by e.g. iLeach fc PetrosianI (119831 ). With the observational 
capabilities of instruments such as Yohkoh and RHESSI it became possible to directly com- 
pare the models with observations. iPetrosian fc Donaghyl (119991 ) investigate the conditions 
for formation of loop-top sources and c ompare predicted time- profiles of the X-ray emission 
with time-profiles observed by Yohkoh. iMinoshima et al.l (120081 ) use a trap-plus-prec ipitation 
model to explain RHESSI observations. While those studies all focused on spectra, iFletcher 
( 119961 ) investigated the height of X-ray sources obtained from numerical simulations and 
compared them with Yohkoh observations, finding that they are consistent with partial elec- 
tron trapping in a magnetic trap. However, there are no quantitative studies of the HXR 
source sizes. 

In this paper, we focus on the size of HXR sources produced by non-thermal electrons 
{E > 30 keV) in Coulomb collisions with the ambient plasma. Within the assumption of 
a single monolithic loop we analyze a range of processes which could increase the vertical 
extent (along the magnetic field lines) of the sources and explain the observations. Thus we 
consider: i) density variations in a single loop; ii) the role of the initial pitch angle distribution 
iii) the effect of pitch-angle scattering; iv) magnetic mirroring, and v) instrumental effects 
related to RHESSI observations. 



2. Transport of energetic electrons in X-ray footpoints 

The evolution of an initial electron flux distribution Fo{Eq, z) as a function of distance z 
alon g the magnetic field lines is described by the time independent Fokker-Planck equation 
(e.g. iHolman et al.ll201ll . as a recent review), including magnetic mirroring and coUisional 



pitch-angle scattering and energy loss: 



dF 

dz 



1 - /i^ d(ln B) OF Kn{z) dF n{z)K d 



2/i dz dfi fiE dE 2E^fx d^i 



OF 



Kn{z) 
' fiE^ 



where n{z) is the plasma density at distance z, K = 27re^A, with A the Coulomb logarithm, 
B is the magnetic field strength, and fi is the cosine of the electron pitch angle relative to the 
magnetic field direction. In the case of purely coUisional energy loss in a uniform magnetic 
field and an injected electron power-law fiux distribution Fq{Eq) ~ E^^ this leads to: 



F{E, z) = Fo[Eo{E, z)]E.jE^ + 2KN{z 



-(-5-1) 



(2) 
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where Eq{E,z) = {E'^ + 2KN{z)Y^'^ and N{z) = J n{z)dz is the column depth. This we will 
refer to as "simple thick-target" in this paper. Equation [1] can also be sol ved analytically by 



neglecting collisions but inclu ding the term for the magnetic field (e.g. iLeach fc Petrosian 
198ll : IZweibel fc Haberlll983l ). However, this analytical solution is only applicable as long 
as the electron pitch angle cosine is larger than about 0.7 because of a factor ~ l//i that is 
introduced in Eq. [2l Thus, if one wants to include all physical effects one has to use test- 
parti cle simulations. Equation [1] can be writt en as a set of stochastic differential equations 



(e.g. iKovalev fc Korolevlll981t iFletcherl 119961 ). Since the typical integration time necessary 
for reliable imaging of footpoints with RHESSI is of the order of 30 - 60 seconds, which 
is very large compared to the electron loop transit time and collisional loss time in the 
chromosphere, a time-independent treatment is a good approximation. Therefore we consider 
time-independent equations along the particle path: 



dz 

ds 
dE 

ds 

dfi 

ds 



dlnB fiKn{z) 



2dz 



E^ 



+ 



Kn{z) 



nl/2 



(3) 
(4) 

(5) 



where z is the distance from the point of injection along the magnetic field lines, s is the 
path of the electron, and W{s) denotes a standard Wiener process. Equations [3H5] are solved 
using a numerical scheme: 



Zj+l 

Ej+i 



Zj + fijAs 

[E] -2Kn{zj)AsY^^ 



(6) 
(7) 



dlnB 



dz 



As 



fijKn{zj) 



As 



,,Kniz^ 



1/2 



where ^ is a random variable taken from the normal distribution p(^) = l/y/27iexp (— ,^^/2) 
for each step As. The schem e proved to be reliable in modelling the effect of collisional 
scattering on pitch angle (see iBail Il982l : [MacKinnon &: Craigl Il99ll ). Equations ([3IH1) are 
solved using a power-law electron distribution 



FoiEo,h = h 



loop) 



E~. 



(9) 



injected at the top of the loop with high energy cutoff energy at 500 keV and injection 
height hioop = 12.5 Mm. Equation [2] was used to test the simulations. Both numerical and 
analytical solutions result in the number of electrons as a function of distance z from the 
injection point, or, equivalently, as a function of height h above the photosphere F{E, h) 
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where h = hioop — z. Finally, calculating the corresponding X-ray bremsstrahlung intensity 
/(e, h), and determining the X-ray source vertical profiles, this is compared to observations 
made with RHESSI. 



2.1. Definition of size and position 



In the observations of limb events presented in lKontar et al.l (I2008l . l2010l ) and lBattaglia fc Kontar 



(I2OIII ). we found the position and size by forward fitting X-ray visibilities (2D spatial Fourier 
components) with a circular or elliptical Gaussian source model. The size of a source is 
thereby defined as the full-width at half maximum (FWHM) of the fitted Gaussian. 



27r(Tr(Ti 



exp 



(x-Xo(e))2 (y-yo(e))^ 



2^2 



(10) 



where 2 a/2 \ia2ax and 2 a/2 \n2ay are FWHMs of an elliptical Gaussian source in the x and 
y direction respectively. For HXR sources observed at the solar limb, the size along the 
radial direction represents the vertical FWHM size of the source, while the perpendicular 
size (along the solar limb) is equivalent to the size of a HXR footpoint parallel to the solar 
surface. The simulations performed here provide the electron fiux as a function of height 
F{E, h) (one dimension) and hence the vertical extent of HXR sources. The X-ray fiux 
/(e, h) per unit distance is then given as: 



I{e,h) 



n{h)A{h) 



F{E,h)a{e,E)dE 



where A{h) is the area of the magnetic fiux tube at height h, R is the distance Sun-Earth 
and a{e, E) the bremsstrahlung cross-section. 



We define the height of a source as the first moment of the X- 
FWHM as the second moment. Only emission larger than 10% of the 
to com pute the mornents, to emulate the fact that RHESSI images 
range (IHurford et al.l |2002| ) . Figure [1] illustrates the photon fiux as 
two photon energies in the analytical solution of the simple thick-targ 
side of the figure shows the observed FWHM in the event analyzed 
compared with the expected FWHM in a simple thick-target model 
solution of Eqs. PITT]) . 



ray fiux profile and the 
maximum fiux was used 
have a limited dynamic 
a function of height for 
;et case. The righ t -hand 



m 



Kontar et al.l fl2010f ) 



given by the analytical 
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Fig. 1. — Left: Normalized photon flux as a function of height above the photosphere at 30 
keV and 70 keV (red), using a density profile of n{h) = 10^° + 1.16 x 10^'' x exp{—h/ho) 
(comp. Eq. [T2|) with a scale-height ho=130 km. The vertical dashed lines give the first 
moment (maximum position) of the emission above 10 % (horizontal dashed line). Right: 
Observed FWHM in the event of January 6th 2004 (after Ixontar et al.l[201ol ). The black 
line is the expected source FWHM in a simple thick target, using Eqs. ElandfTTl 
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3. Analytical and numerical results 

We first investigate tlie effect of the cliromosptieric density function on the resulting 
source size in the simple thick target (analytically), using Eq. |2l Then, the influence of the 
initial pitch-angle distribution, coUisional pitch-angle scattering and magnetic mirroring will 
be explored using test-particle simulations (Eq. EJ-E]). Finally, in Section 1X71 we will discuss 
in more detail how the sources would be observed by RHESSl and what results we would 
expect using visibility forward fitting. 



3.1. Chromospheric density structure 



The deposition of large amounts of energy into the chromosphere by energetic elec- 
trons could lead to a number of processes, including heating and expansion of the chro- 
mosphere. This process of chromospheric evaporation generally leads to redistribution of 
plasr aa density in the flaring atmosphere, increa sing the density of plasma in the flaring 
loop (lHirayamalll974j : lAntiochos fc Sturrockl 119781 ). At the same time, the density structure 
in a flaring loop can strongly affect the source size even in the case of purely coUisional 



energy loss. In previous work, a single scale-height exponenti al densitv (IKontar et al. 



2010l : iBattaglia fc Kontarll201l(). multiple scale-heigh t density ( jSaint-Hilaire et al.ll2010[ ) or a 



2008 



power-law density function ( lAschwanden et al.ll2002l ) have been investigated. Although the 
positions of HXR sources are in agreement with the single scale-height exponential density 
profiles with scale heights of 130-200 km, the predicted v ertical source s izes a re up to a factor 
of 4 smaller than observed (IBattaglia fc Kontaiil201ll ) (IKontar et al.l ( 120101 ) suggested that 
a multi-threaded density structure with vertical strands of different density could increase 
the source size significantly). The exponential density function 



n{hyP = ni + riphot x exp{-h/ho) 



:i2) 



is used by IBattaglia fc Kontan (120111 ) where n j, is the (constan t ) loop density, Uphot = 1-16 x 
10^'' cm~^ the photospheric density following IVernazza et al.l ( 19811 ). and Hq is the density 
scale height. Figure H] illustrates such an exponential density model with scale-height /io=130 
km, and also a density model of the shape of a K-function 



n. 



phot 



X 



h 



(13) 



for a K, factor of 10 and two different scale heights ho = 130 km and ho = 290 km (dark blue 
and light blue curves in Fig. |2]). In addition, densities that combine the exponential function 
near the photosphere with either a quadratic function or a /t-function at larger heights were 
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used, where the exponential plus quadratic function was given as n{h) = n{hY^P + 1.16 x 
10°'(/Zmax — where h^ax is the loop height and a = [10, 11, 12] (yellow, dark red, and 
red curves in Fig. [2]). The exponential plus K-function was given as n{h) = n{hY^P + n{hY 
(green curve in Fig. |2]). These latter models account for the fact that the high density in the 
lower chromosp here, below about 1 M m, will not be affected significantly by processes such 
as evaporation (iMachado et al.lll980l ). while the loop density may be considerably different 
from quiet Sun density models. As Fig. [2] demonstrates, a notable effect on the FWHM 
is only found if the loop density reaches extreme values of more than 10^^ cm^^. Such 
coronal densities a re rather non-typically high although are observed in some flares (e.g. 
Veronig et al.ll2006l ). Such high densities will lead to energetic electrons of energy ~ 20 keV 
and even higher being coUisionally stopped in the coronal part of the loop before reaching 
the chromospheric footpoints. 

All these density functions result in an effective increase of the loop density compared 
with the single scale-height exponential and the resulting HXR source FWHM is larger by 
up to a factor of 4 than in the case of an exponential density (Fig. [2]). However, the location 
of the peak of the emission is also found to be a factor of 4 higher, and the bulk of the 
emission below 40 keV comes from the top of the coronal loop, therefore no footpoints would 
be ob served below 40 keV, which is not the case in the observations ( iBattaglia fc Kontar 
20111 ). where footpoints are observed at energies as low as 20 keV. 






Height [Mm] 



Energy [keV] 



Energy [keV] 



Fig. 2. — Effect of different density models on the height and FWHM of HXR sources. 
Left: Density models. The black curve is an exponential density with a single scale- height 
of 130 km. Dark blue and light blue are ^-functions with different scale heights. Green is 
an exponential density combined with a /t-function. Yellow, red and dark red represent an 
exponential density combined with a ~ /i^ function (see text for details). Middle: Height of 
peak emission as a function of energy. Right: FWHM as a function of energy. 
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Due to the relatively modest increase of the chromospheric temperature at heights below 
1 Mm, we assumed a neutral atmosphere. However, the plasma in the transition region and 
lower corona will be partially or completely ionized with a change in the ionization fraction 
at some height in the chromosphere. The ionization state of the medium affects the Coulomb 
logarithm, so that InAejy ~ 7.1 in neutral media and In Age — 20 in fully ionized plasma, 
where Aee and KeH are the electron-electron and electro n-hydrogen C oulomb logarithms. 



The effect of this on HXR s pectra has be e n disc ussed by iBrownl (119731 ) and in the context 



of RHESSI spectroscopy by iKontar et al.l (120021 ). To investigate the influence of ionization 



change on the height and the FWHM of HXR sources, we introduced an effective Coulomb 
logarithm 

A = (Aee - Ae^) [x{h) + J^^^^^ , (14) 

The ionization fraction x{h) was assumed to be a step function with x{h) = 1 for h > 0.9 
Mm and x{h) = for h < 0.9 Mm. The effect of non- uniform density is illustrated in Fig. 
[3] along with a comparison with the case of constant Coulomb logarithms. The case of a 
more realistic atmosphere with partial ionization would lead to results within the extremes 
illustrated in Fig. |3] and will not change the source size noticeably. Thus, the use of an 
idealized ionization structure is justified. 

The simulations show that although different chromospheric and loop density models 
can increase the source size by up to a factor of 4 to fit the observations, this will also 
change the height of the sources, contrary to the observations. It has to be noted that in 
order to produce a noticeable change of the HXR source size, the density structure of the 
whole atmosphere, including the transition region and corona needs to be changed by a few 




Energy [keV] Energy [keV] Energy [keV] 



Fig. 3. — Left: Total photon spectra for completely neutral target (black), completely ionized 
target (blue) and ionization step-change (red). Middle and right: Height and FWHM of HXR 
sources as a function of energy in the three cases. 



orders of magnit ude. This seems in contradiction with both theoretical models of the flaring 
atmosphere (e.g. lAUred et al.ll2005l ) and observations. 



3.2. Pitch-angle distribution 

The height and the width of a source is also likely to depend on the initial pitch-angle 
distribution of the injected electron beam. In the simple thick-target model, injection and 
propagation of energetic electrons is assumed parallel to the magnetic field (i.e. cosine of 
initial pitch angle /io = 1), but in the other extreme case of injection perpendicular to the 
field (/io = 0) or strong trapping, the electrons would lose energy near the acceleration 
region, the height would be constant as a function of energy, and the FWHM would depend 
on the extent of the acceleration region. In an intermediate situation < /io < 1, energetic 
electrons are expected to lose energy at different heights depending on their initial pitch- 
angle distribution. The effect of an initial pitch-angle distribution is investigated including 
collisional energy loss, but no change in the initial pitch-angle distribution is assumed, i.e. 
Eq. |5] is dfi/ds = 0. In this and the subsequent sections, an exponential density profile 
n{h) = rii+nphot exp(— /i/Zio) with a scale height of 144 km and a constant Coulomb logarithm 
of A = 7.1 is used. This density profile is consistent with the source heights measured with 
RHESSI and the theoretical modeling of the low atmosphere at heights h <1 Mm. 

Figure S] illustrates how the initial pitch-angle distribution affects the height and the 
FWHM. Three different cases are presented: /io = 1 (as a reference), 0.9 < /io < 1 uniform 
(strongly beamed) and 0.1 < /io < 1 uniform (broad distribution). The broad distribution 
leads, as expected, to a larger source size as a function of energy, as well as a larger FWHM. 
However, the maximum change in both height and FWHM is about 10% in the case of the 
broad distribution. Therefore, this effect alone cannot account lead to the observed vertical 
sizes of HXR sources. 



3.3. Collisional pitch-angle scattering 

In collisional interactions with the ambient plasma, electrons do not only lose energy, 
but are pitch-angle scattered with a similar rate. Therefore, their pitch-angle distribution 
changes as the particles propagate downwards towards dense regions of the atmosphere. 
Figure |5] compares the standard case (collisional energy loss only), with the outcome of a 



- 11 - 



situation with collisional pitch-angle scattering included, i.e. Eq. [5] becomes 



ds 



fiKn{z) 



+ 



Kn{z) 



1/2 



w. 



(15) 



In this case, the maximum emission is located up to 20 % higher than in the no-scattering 
case. However, the effect on the FWHM is small, at the level of about 5%. As expected, the 
effect is qualitatively similar to the injection of a broad initial pitch-angle distribution. 



3.4. Magnetic mirroring 

A converging magnetic field at the loop footpoints also changes the pitch-angle distri- 
bution of energetic electrons, causing electrons with large pitch angle (small cos 9) to mirror 
upwards from the footpoints before they are collisionally stopped. This might further con- 
tribute to an increase in source size. In this Section we consider collisional energy loss, 
and pitch-angle change due to magnetic field convergence, but not collisional pitch-angle 
change, i.e. Eq. [5] becomes dfx/ds = —(1 — A'-^)^^^- The magnetic field strength is mod- 
eled as B{h) = Bq + i?itanh( — (/i — huMkM), which adequately represents a converging 



magnetic field in the chromosphere (IFedun et al.ll201ll ). This model gives a magnetic field 



B{h » hnj) Bq — Bi at coronal heights and B{h = 0) — t- i?o + -Bitanh(l) at the pho- 
tospheric level. The increase in field strength relative to the ambient density is illustrated 
in Fig. ini The field convergence (and the electron pitch angle) defines the depth of the 
mirroring point. If the magnetic field converged higher up in the loop, the sources would 




Fig. 4. — Normalized photon flux as a function of height at 30 keV and 70 keV (solid 
and dashed lines, left), height of maximum emission and FWHM as a function of photon 
energy (middle and left), for three cases of initial pitch-angle distribution: /io=l (black), 
0.9 < /io < 1 (red), 0.1 < /iq < 1 (blue). 
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Fig. 5. — Height of maximum emission and FWHM as a function of photon energy (middle 
and left), for two cases of initial pitch-angle distribution: 0.9 < yUo < 1 (red), 0.1 < fio < I 
(blue), and including collisional pitch-angle scattering. The black line indicates the result 
for collisional energy loss only. 

be observed higher up. In the extreme case, one could simulate a coronal source caused by 
electron trapping. 




12 3 4 5 100 100 

Height [Mm] Energy [keV] Energy [keV] 



Fig. 6. — Left: Density and magnetic field strength as a function of height. Height of 
maximum emission and FWHM as a function of photon energy (middle and right), for two 
cases of initial pitch-angle distribution: 0.9 < fJ^o < I (red), 0.1 < fJ^o < I (blue), and 
including magnetic mirroring. The black line indicates the result for collisional energy loss 
only. B{h) = Bq + Bitanh{-{h - hM)/hM) with Hm = 1 Mm, Bq = 600 Gauss, Bi = 500 
Gauss. 
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In the case presented in Fig. [6l the source height is increased by a factor of ~ 1.6, while 
the FWHM increases by a factor of 1.7 - 1.3, depending on energy. Although the size increase 
is larger than in the case of coUisional scattering, it is still not strong enough to explain the 
FWHM observations. In addition, a nearly isotropic initial distribution of electrons (Fig. E]) 
leads to a larger source FWHM at higher energies, contrary to X-ray observations. 

3.5. Magnetic mirroring and collisional pitch-angle scattering 

Finally, the effects described in the above two sections are combined and the full Eqs. 
OlSlare solved numerically (Fig. [7]). Since the effect of pitch-angle scattering itself is small 
compared to the effects of magnetic mirroring, this case is dominated by the effect of the 
magnetic field and the result very similar to the case of magnetic mirroring (Fig. [7]). 




Energy [keV] Energy [keV] 

Fig. 7. — Height of maximum emission and FWHM as a function of photon energy (middle 
and left), for three cases of initial pitch-angle distribution: /io=l (black), 0.9 < /io < 1 (red), 
0.1 < /io < 1 (blue), and including both collisional pitch-angle scattering and magnetic 
mirroring. 



3.6. Constructing an electron distribution as a function of height 

The processes described in Sections 13.31 - 13.51 all influence the electron distribution as a 
function of height, e.g. pitch-angle scattering causes electrons to be stopped higher up in 
the loop. However, as illustrated in Fig. E] - [71 this has a rather small effect on the X-ray 
flux profile and thus the source FWHM. The FWHM of an X-ray source is proportional 
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to the product of electron flux density and plasma density I{e,h) ~ F{E, h)n{h). In the 
extreme case of F{E,h) ~ l/n{h), the resulting I{e,h) is independent of height and this 
constant value could extend vertically over all h where F{E,h) ~ l/n{h). We can therefore 
ask how F{E, h) should look in order to make the product of F{E, h)n{h) independent of h, 
hence increasing the size of the X-ray source. Starting with F{E, h), as found in the simple 
thick target case, we modified the shape of F{E, h) for every energy, so that the slope of the 
curve was close to l/n{h), as shown in Fig. [8] (top left), then we computed the height of the 
maximum emission and FWHM. As F{h) approaches l/n{h), the FWHM of the resulting 
X-ray flux increases up to 5 arcsec. 



10" 




. — 1 0''5/density _ 



1 ? 3 

Height [Mm] 





Energy [keV 



Fig. 8. — Top left: Electron flux as a function of height for 24 keV electrons. The solid 
line represents the simple thick target case. The dashed and dotted lines are "artificial" 
distributions, constructed so that the slope nears that of l/n{h). The blue dash-dotted line 
illustrates l/n{h). Top right: normalized X-ray flux. Bottom left: height and FWHM (red 
curves) for the three different electron distributions where the line styles correspond to the 
line styles in the top left panel. Bottom right: Injected electron spectrum (black line) and 
electron spectrum at a height of 0.8 Mm. 
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However, the height of the resulting HXR maximum as a function of energy is constant, 
as is the FWHM. Further, the electron spectrum at low heights is completely different from a 
thick target spectrum (Fig. El bottom right). Most importantly, such an electron distribution 
would have to be extremely "fine tuned" to the ambient density. 



3.7. Instrumental effects and method 

Comparing the modeling of the main electron transport effects with observations and 
finding physical explanations for the observed source sizes, we assumed that the observed 
difference is not entirely due to instrumental effects. This was based on there being no 
modulation in the finest RHESSI grids (Grid 1 has a spatial res olution of ~ 2.3 arcsec an d 



grid 2 has spatial resolution ~ 3.9 arcsec) in the observed events (IBattaglia fc Kontaiil201ll ). 
indicating that the source dimensions must be of the order of > 4 arcsec. Here we perform a 
more quantitative study of the instrumental effects, using the simulation software developed 
by Richard Schwartz (private communication). The simulation software uses an arbitrary 2 



dimensional map as input and creates a corresponding calibrated event list (ISchwartz et al. 



20021 ) . The standard RHESSI imaging algorithms are then used to construct the image. 
We used several source models (circular Gaussian, elliptical Gaussian) as the initial map 
and forward-fitted the visibilities from the corresponding calibrated eventlist to compare the 
fitted FWHM with that of the original map. We find that circular Gaussians are correctly 
recovered within the uncertainties, e.g. a circular Gaussian with FWHM 1 arcsec is fitted 
with a FWHM of 0.93 ±0.1. In the case of an elliptical Gaussian, the fitted major and minor 
axes tend to be larger than the original ellipse, e.g. an input elliptical source with major and 
minor axes of 3 and 1 arcsec respectively is fitted with 3.1 and 1.5 arcsec. Thus the fitted 
minor axis is 50% larger than the input minor axis. 

This last example shows that, while there are instrumental effects, especially in the case 
of elliptical sources, these effects cannot entirely account for the observed source sizes. 

The FWHM found in the previous sections are all in terms of the second moment of the 
1-dimensional height distribution of X-rays I{e,h). However, the limited dynamic range of 
RHESSI influences the accuracy of the measured moments. X-ray flux in RHESSI images 
which is around or less than 10% of the brightest part of the image is dominated by an error 
from the brightest source. However, using the simulations we can address the question of 
how such a source would be observed with RHESSI, and how the second moment found in 
the simulations relates to an observed RHESSI image. As input we used the shape of the 
photon-fiux as a function of height found in case of a density n{h) = n{hY^^ + n{hY (green 
line in Fig. [2]) with 6 arcsec width. This map was used as input for the simulation software 



- 16 - 



and Clean and Pixon images (using software defaults and grids 1-8) were reconstructed from 
the calibrated event lists. Figure [9] displays the profile of the photon flux as a function of 
height, the input map, and the Clean and Pixon maps. The resulting fitted FWHM are 4.9" 
and 2.7" for the major and minor axes. The second moment of the 1-dimensional X-ray flux 
distribution are is 2.8" and becomes 2.1" if computed only for the emission exceeding the 10% 
level. Therefore, the second moment of the full distribution overestimates the size, while the 
moment of the flux > 10% underestimates the size, compared to visibility forward fitting. 
It has to be added that the source model was very simple and there is no unmodulated 
background added to calibrated event lists. 



4. Summary and conclusions 



In the simple collisional thick-target model, both the height of the maximum HXR 
emission and the vertical HXR source sizes are determined by the density scale height. 
RHESSI observations suggest that the HXR source positions can be well fitted with a 
single exponential scale-height density model, assuming a simpl e collisional thick target . 



This results in scale heig hts between about 130 km and 200 km (IKontar et al.l 



Vernazza et al 



Battaglia fc Kont ad 20111 ). consistent with chromospheric models (e.g. iMachado et al. 



2008 



2010 



1980 



198ll ). However, the observed HXR sizes are about 4 times larger than ex- 



pected from the simple collisional transport model. Here we have quantitatively investigated 
how the density profile, collisional pitch-angle scattering, magnetic mirroring, as well as 
instrumental effects affect the source sizes. 



In lBattaglia fc Kontarl (120111 ) we showed that projection effects and source motion over 
the RHESSI image time interval cannot account for the observed source sizes. In the present 
work, applying RHESSI visibility forward fitting on simulated HXR source maps we demon- 
strate that the source size cannot be due to instrumental effects, alone. This leaves the 
physical effects of magnetic mirroring and collisional pitch-angle scattering which we inves- 
tigate by solving the Fokker-Planck equation both numerically and analytically. While pitch 
angle and magnetic mirroring effects change the electron flux distribution, these effects tend 
to increase the FWHM of the X-ray source profile by only up to a factor 1.5, which is not 
enough to explain the observations. The dominating factor that determines the X-ray source 
size is the atmospheric density structure. In the case of an exponential density model with 
a single scale height and a constant coronal loop density of around lO^'' cm~^, the X-ray 
emission will originate predominantly from the region of highest density. Thus, even though 
the effects described above alter the electron distribution as a function of height, emission 
by electrons higher up in the loop will always be faint compared to the emission from the 
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Fig. 9. — Top left: X-ray (black) profile for simulations with high loop density (blue). Top 
right: Simulated map of a loop with ~ 6 arcsec width and the height profile of the simulated 
X-ray profile. Bottom: RHESSI CLEAN and Pixon image representation of the simulated 
map. The contours represent the 50%, 70%, and 90% levels of the elliptical Gaussian found 
with visibility forward fitting. 

denser chromosphere. Source sizes of around 4 arcseconds can only be achieved by unlikely 
loop densities of the order of 10^'^ cm~^. Such high densities will also cause the HXR sources 
to appear at larger heights, well above typical chromospheric heights. Thus, within the 
traditional thick target model, the only plausible explanation for the observed HXR source 
sizes remains a multi-threaded density structure. 
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